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ABSTRACT 

We investigate the evolution of coronal loop emission in the context of the 
coronal magnetic field topology. New modeling techniques allow us to investi- 
gate the magnetic field structure and energy release in active regions. Using these 
models and high resolution multi- wavelength coronal observations from the Tran- 
sition Region and Coronal Explorer (TRACE) and the X-ray Telescope (XRT) on 
Hinode, we are able to establish a relationship between the light curves of coronal 
loops and their associated magnetic topologies for NOAA Active Region 10963. 
We examine loops that show both transient and steady emission, and we find 
that loops that show many transient brightenings are located in domains associ- 
ated with a high number of separators. This topology provides an environment 
for continual impulsive heating events through magnetic reconnection at the sep- 
arators. A loop with relatively constant X-ray and EUV emission, on the other 
hand, is located in domains that are not associated with separators. This result 
implies that larger-scale magnetic field reconnections are not involved in heating 
plasma in these regions, and the heating in these loops must come from another 
mechanism, such as small-scale reconnections (i.e., nanoflares) or wave heating. 
Additionally, we find that loops that undergo repeated transient brightenings are 
associated with separators that have enhanced free energy. In contrast, we find 
one case of an isolated transient brightening that seems to be associated with 
separators with a smaller free energy. 

Subject headings: Sun: corona — Sun: magnetic topology — Sun: UV radiation 
— Sun: X-rays 
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Introduction 



Observational studies of the character and evolution of coronal structures are becoming 
more refined as observational capabilities improve and accompanying models by which to 
interpret the results are developed. In studies of the evolution of coronal emission, obser- 
vations at coronal-emission wavelengths show a variety of structures that evolve in distinct 
ways. It is useful to examine these observations in the context of the magnetic field. Coro- 
nal loops observed in soft X-ray and EUV have been stud i ed since the Skylab observations 
in the 1970s (e.g., IVaiana et al.lll973l ; lYoshida et al.lll995t iBerger et al.lll9991 ). In previous 
studies, the spatial and temporal correlations between soft X-ray observations and EUV ob- 
servations of coronal loops have been investigated in o rder to study coronal loop ev olution 
dSchmieder et al.lbood ; iNittalboOol : iNagata et aDl2003h . Iwinebarger fc Warren! tood ) found 
that structures observed in the X-rays and then subsequently in the EUV are due to emission 
from heated coronal loops cooling through the passband of each instrument. 

Recent studies using the Hinode/EXJV imaging spectrometer (EIS) have presented ob- 



servations that show a brightening in EUV following a brightening in X-rays (j Warren et al. 



20071 ; lUgarte-Urra et al.ll2009t IWarren et al.ll2009l . 120101 ) . A multi-thread hydrodynamic sim- 
ulation has reproduced several coronal loop characteristics, including high electron density 
and a long observed life time, for the coronal loops near 1 M K, but it has difficu lty reproduc- 
ing the high temperature loop emission observed by XRT ( IWarren et al.l 120091 ). It has also 
been shown that sim ple impulsive and quasi-s teady heating does not reproduce the observed 
loop characteristics ( lUgarte-Urra et al.l 120061 ). 



Sever al studies have shown evidence for magnetic field changes during the coronal loop 



evolution (IWinebarger & Warrenl l2005t IWarren et al.l 120071 ). Total magnetic flux in active 
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the field strength and its variation at loop boundaries. One recent study suggests that the 
changes in the topology of the magnetic field may be related wit h the evolution of the high 



temperature coronal loop in X-ray emission ( IWarren et al.l 120091 ). 



see 



Starting from m agnetogram observations, Magnetic Charge Topology (MCT) models 



Longcopell2005l . for an overview) provide a framework to analyze the coronal topology, 
including the magnetic field connectivity and separator field lines. In this class of models, flux 
concentrations on the boundary are represented by point charges or sources. The magnetic 
field lines established by the MCT model form domains, de fined as continuous yolume s 
containing field lines that connect the same pairs of sources ( jLongcope fc Klapperl 120021 ). 
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Thus, the MCT model can determine which domains are along a given line of sight in a 
potential (current-free) extrapolation of the coronal magnetic field. 

Domains are bounded by separatrix surfaces, whose intersect ions form se p arato r field 
lines, and thus are the location of reconnection in MCT models. iPriest et al.l ( 120051 ) have 
proposed that slow photospheric motions can generate electric current near separa tors and 
separatrices which can be released as a heating energy. Recently, iNoglik et al.l (120091 ) showed 
that a potential magnetic field extrapolation represents the structure of the hotter X-ray 
loops and the larger cool loops seen in 171 A i mages appear to follow the separatrix surfaces 
determined by the MCT analysis. In addition, iPlowman et al.l (120091 ) show that separators 
explain the observed uniform width of coronal loops better than randomly distributed po- 
tential field lines. In this paper we use the MCT model to determine the location of flux 
domains and separators. By construction, this model disallows reconnection, but can iden- 
tify the location of the separators where reconnection is most likely to occur. The details of 
how reconnection across the separators heats coronal loops in a flux system requires a more 
detailed time dependent MHD simulation that is beyond the scope of this paper. 



The Minimum Current Corona (MCC) model (ILongcopd 119961 ) has been developed in 
order to determine the quasi-static evolution of magnetic fields starting with the potential 
fields of the MCT model and introducing currents along separators in response to footpoint 
motions. The MCC model represents the corona with the minimum permissible current after 
the evolution usi ng a series of qu asi-static steps and subject to the constraint that there is 
no reconnection (lLongcopdll996l ). The resulting increase in en ergy above the potential field 
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energy (fre e energy) can be released in the form of heating (ILongcopd Il99 
2000l . 120051 ). This model has been applied to several flare events to investigate the energy 
source of the flares (ILongcope et al.l 120071 ; Ides Jardins et al.l 120091 ). 



In this paper, we use the MCC model to calculate the energy build-up from the photo- 
spheric magnetic field evolution, presenting this output in the context of mult i- wavelength 
observations of the temporal evolution of the coronal emission. The MCT model is used to 
establish potential magnetic field connectivity, determine separators and domains along the 
line of sight at various coronal heights. The MCC model is used to determine coronal cur- 
rents and quantify the energy buildup. The results of the MCC analysis are compared to the 
evolution of the coronal loop brightness, comparing selected areas of both steady emission 
and transient emission. 

In §2, we describe the observational data employed in this analysis. In §3, we explain 
briefly the MCT and MCC models. We also describe the magnetic topology of AR 10963 as 
characterized by the MCT model and the free energy calculated by the MCC model in §3. 
In §4, we discuss the observed coronal loop brightness evolution as compared to the results 
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of the MCT and MCC models. In §5, we present our conclusions. 



Observational Data 



The target of this study is NOAA Active Region 10963 (AR 10963), which transited 
the visible solar disk from 2007 July 8-23. This region is a very small, exceptionally 
quiet numbered active region, producing a few B- and small C-class flares (peak flux in 
the Geostationary Operational Environmental Satellite (GOES) 1-8 A soft X-ray band of 
~ 10~ 7 , 10- 6 Wm 2 , respectively). The region was event-quiet above the B-class threshold 
from July 11 00:05 UT until July 18 13:50 UT. The day of interest for this study was during 
the quiet time, focusing on 2007 July 14. 



The X-ray telescope (X RT, iGolub et al.ll2007l ) on board the joint NASA/ JAXA mission 
Hinode (IKosugi et al.ll2007l ) observes the solar corona at high (~ 1" pixels) spatial resolution 
in multiple bandpasses. For this study, we use a sequence of images from the Al/Poly filter, 
which samples plasma at > 2MK, with a 20 s cadence from 18:27 UT on July 14 to 01:47 UT 
on July 15. The XRT observed AR 10963 with deep exposures in the C/poly (> 3MK) filter 
to study faint structures until 18:15 UT on July 14. These earlier observations do not allow 
for study of the changes in coronal brightness since coronal loop images were saturated. 



The Transition Region and Coronal Explorer ( TRACE, lHandy et al.lll999l ) observations 
are at a 1 minute cadence in the 171 A band, which samples the coronal plasma at ~ 1 MK 
with ~0.5" pixels. 

The line-of-sight component of the pixel-averaged photospheric magnetic field was ob- 
tained from the Mich aelson Doppler Inter ferometer (MDI) aboard the Solar and Heliospheric 
Observatory (SoHO, IScherrer et al.lll995[ ). For this study, a roughly 24 hour time sequence 
of the 96-minute cadence, ~ 2" data was used, from 2007 July 14 01:39 UT to 2007 July 15 
00:00 UT. 



2.1. Co-alignment among XRT, TRACE, and MDI 

Co-alignment between the datasets arising from three different instruments on three 
different spacecraft was crucial for this study, and was performed in the following manner. 

We used the MDI level 1.8 data and full-sun, level 1 data from the Extreme-ultraviolet 
Imaging Telescope (EIT) on SoHO (calibrated by eitjprep.pro in the SolarSoft suite of analy- 
sis tools). We first supposed that an internal co-alignment between MDI and EIT is correct, 
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and the EIT images then provided the common coalignment platform for both the XRT 
and TRACE data. Active region 10963 was observed by XRT and TRACE for ~6 hours 
with Al/Poly (1.0286" pixels) and 171 A (0.5" pixels), respectively. The EIT (2.63" pixels) 
observed the full sun every 6 hours at 284 A and 171 A. 

The XRT observed the full solar disk with Ti/poly every 20 minutes as a context image. 
X-ray bright points on the full Sun at 19:07:15 UT on July 14 were aligned with EIT 284 A 
at 19:06:07 UT on July 14. Due to the similar response to high temperature plasma, the EIT 
284 A and the X RT Ti/poly observat ions were used to co-align the satellites. The spacecraft 



jitter on Hinode (jShimizu et al.ll2007l ) was corrected with the reference image at 19:07:15 UT. 



The TRACE 171 A observation at July 14 18:59:26 UT was aligned with the EIT 171 A 
observation at July 14 19:00:13 UT. The aligned image was used as a reference to align for 
later time observations for ~ 6 hours (using trace_cube_pointing.pro in SolarSoft) to account 
for jumps in the TRACE pointing as the active region was tracked across the disk. 



2.2. Area selection for light curves 

We select four areas in the XRT and TRACE data with which to investigate coronal 
brightness changes. These areas are presented as four boxes in Figure [U Inner and outer 
boxes correspond to approximately 5" x 5" and 10" x 10" areas, respectively. The Boxes 1, 
2, and 4 were selected based on the coronal loops in the XRT observation, and Box 3 was 
selected based on a small loop visible in the TRACE data at 22:24 UT on July 14. 

Light curves for each box are presented in Figure [21 Box 1 shows a brightening in the 
TRACE observation about an hour and a half after a decrease in X-ray brightness. Box 2 
shows transient brightenings in XRT, however the TRA CE light curve shows only very small 
changes in brightness. Box 3 shows transient brightenings in both XRT and TRACE. As seen 
in the animation version of Figure [U the light curves of Boxes 2 and 3 include the brightness 
changes for several loops. Box 4 shows steady emission in both XRT and TRACE, at a 
lower level in XRT than any of the other boxes. The light curves in Figure [2] show similar 
behaviors for the two box sizes, which verifies that the light curves do not depend on a small 
change in the location of the box. 

Figure [T] also shows a loop (LI) that contributes to the light curves in Box 1 on the 
XRT and the TRACE observations at 22:24 UT and 23:48 UT, respectively, when the bright- 
enings are seen in each instrument. The ends of the loops seen by the XRT observation are 
represented with Fl. These locations are used to estimate the height of the loop in §4.1 (see 
also Table H] in Appendix). The location of the loop for the Box 1 is selected based on the 
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loop observed by the TRACE because the loop is better resolved in the EUV observation, 
and then the same location is represented on the XRT observation. 



3. Magnetic Topology and Energy Buildup 

We utilize a Magnetic Charge Topology model, and the Minimum Current Corona model 
to understand the magnetic topology and the magnetic energy changes in the active region 
that we observe with XRT and TRACE. We use the MCT model to define the magnetic 
topology and locate possible sites for magnetic energy reconnection. With the MCC model, 
we quantify the magnetic energy changes. The details of these models and our implementa- 
tion of them for the analysis of AR 10963 are described in the next subsections. 



3.1. Magnetic Charge Topology model 

The MCT model is used to establish the magnetic connectivity by potential magnetic 
field extrapolation into the corona fro m a distribution of point sources representing the pho 



Priest fc Forbes 1989; Lau 1993; Demoulin et al 



tospheric magnetic field distribution flBaum fc Bratenahllll980l; iGorbachey fc Somovl Il988 



1994 



iParnell et al. 



1994 ). As a result of 



the extrapolation, the magnetic field lines form a domain, defined as a volume containing 
field lines that have the same sources at their ends (unless the field lines tend to infinity). 
In the following paragraphs, we explain several topological terms derived from the inferred 
magnetic connectivity that we use in this paper. 

Locations where the magnetic field vanishes are referred to as magnetic null points. In 
three dimensional geometry, the local mag netic structure of a n ull can be described in terms 
of a fan surface and two spine field lines (IParnell et al.lll996l ). The fan surface is a set of 
field lines which radiate out from a null point (positive, B-type null) or radiate into a null 
point (negative, A-type null). The spine field lines are those directed toward the positive 
null point or directed away from the negative null point. 

The fan forms a separatrix surface that divides the corona into volumes of different 
connectivity called domains. Separators are field lines which begin at a positive null and end 
at a negative null, and they are located at the intersection of separatrix surfaces. Magnetic 
reconnection and energy release in a coronal field stressed by photospheric motions are likely 
to occur near separators, although it should be noted that the MCT model does not allow 
for reconnection. 



The flux in a connection (Vfy*) is determined from the number of magnetic field lines 
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initiated from the source (i) i n random directions and followed to their termination at a source 
of the opposite polarity (J) (IBarnes et al.l 120051). For each pair of source s, it is possible to 
have connections through multiple domains (jLongcope &: Klapperll2002l ). In this analysis, 
each different pair of sources defines a single domain. Thus the connectivity flux and the 
domain flux are the same for each pair of sources in this work. 



3.2. Partitioning the magnetogram time-series 

To extrapolate the magnetic field lines into the corona, the MCT model requires a set 
of distinct point sources derived from an observed magnetic field. Therefore, the observed 
magnetic fields must be partitioned in a way that represents the surface magnetic flux con- 
centrations. To determine the flux concentrations from the photospheric magnetic field 
observations, the observed magnetic features are partitioned to preserve important magnetic 
morpholo gy, and the location of each source is determined by the flux weighted center of the 



partition ( IBarnes et al.l 120051 ) . 



We use a time series of 96-minute-cadence MDI magnetograms for the photospheric 
magnetic field observation. First, we select the region of AR 10963 from the full-disk mag- 
netograms for ~ 24 hours, taking into account solar rotation. The selected area retains the 
balance of positive and negative sources within a 4% ratio of the total signed to total un- 
signed magnetic line of sight signal. The selected regions are transformed from the image 
plane (or "plane of the sky") to a plane that is tangent to the Sun's surface at the flux 
weighted center of the selected active region. The "/i-correction" is applied, which assumes 
that all field is radial, and divides the line-of-sight magnetogram by the cosine of the observ- 
ing angle, to approximate the radial field. This process transforms the observed line-of-sight 
pixel-averaged magnetic field component to a tangent plane approximation of the magnetic 
flux, as required by the MCT model. 

The top panel of Figure [3] shows the boundaries which result from partitioning the 
positive and negative sources. We initially determine the flux concentration at the last time 
observed, July 15 00:00 UT, since an earlier study has found that performing the partitioning 
in reverse chro nological order provide s the most consistent partitioning in light of noise and 



solar evolution ( jLongcope et al.ll2007l ). The partitioning parameters used are similar to those 



in Barnes et al. 2005, though they have been adjusted to match this specific active region^. 



1 Specifically, for this dataset we employ a field threshold of 3 times the detection threshold, a smoothing 
depth of 0.5 Mm, a saddle point merging level of 200 G, and a minimum of source flux of 10 4 G Mm 2 (see 
Barnes et al.ll2005l for a description of these parameters). 
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Recently, lLongcope et al.l ( 120091 ) found that there are only slight differences in connectivity 
that result from using different partitioning parameters, so our results are unlikely to be 
sensitive to the choice of partitioning parameters. The bottom panel of Figure E] shows 
domain fluxes among all sources at the same time as the top panel. 

MDI 96-minute data consist of images with integration times of 30 and 300 seconds 
during the 24 hour period of interest, and different integration times could conceivably 
change the inferred flux. To mitigate this effect, the detection level was assigned to be the 
3 a level of a Gaussian fit for the histogram of each MDI observation. This procedure sets an 
appropriate background level for each integration time. No significant systematic difference 
was detected in the resulting fluxes above the 3 a threshold. 

Beginning from the end of the time series, the observed magnetic elements in the mag- 
netograms are a dvected to the earlier MDI observations using Fourier Local Correlation 
Tracking (FLCT, iFisher fc Welschll2008l ). in reverse chronological order for the same reasons 
as mentioned above. The FLCT finds the velocities of features in two successive images by 
computing the cros s-correlation function using standard Fast Fourier Transform techniques 
rtWelschet alJ liooifl These velocities are horizontal velocities for each pixel in the MDI 
observations. The computed velocities are then used to track elements and determine their 
partitioning assignments in subsequent observations. Thus, this procedure preserves the 
partitioning information as consistently as possible through the time series. 

In Figure HI we show a plot of the total unsigned flux for all sources obtained by the time- 
series partitioning, and the time-averaged total flux. It is clear that there is little variation 
during the ~ 24 hours, implying that there is no significant emerging or submerging flux 
during this time period. Therefore, we impose that the flux for each source is also constant 
in time, since the flux changes that do occur are too small to contribute significantly to the 
active region field. We use the time-averaged flux for each source in the MCT model. This 
approach prevents the transfer of flux between neighboring sources in the partitioning, and 
precludes changes in connectivity due to changes in source flux resulting from emerg ence 
or submergence, thus simplifying the usage of the MCC model (lLongcope et al.l 120071 ). as 
described in §3.4. 

The source N4 (near [220", -100"] in Figure |3]) is submerged at 17:36 UT on July 14, and 
only appears in 4 time steps. The average flux of N4 for the four sequential MDI observations 
when it is above threshold is 2.6 x 10 20 Mx. However the average fluxes of other sources are 
in the range of (4.9 — 67) x 10 20 Mx. Therefore, we do not consider the source N4 for this 



2 For the implementation of the FLCT here, we employ time differences betwee n two images of 1. unit o f 
length of a single pixel of 1, and sub- images weighted by Gaussian width of 45, see lFisher fc Welscbl (J2008) . 
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analysis on account of its relatively small flux and short lifetime. This choice helps the 
analysis stay consistent with the assumption that there is no emerging/submerging flux. 

In Figure we demonstrate the evolution of partitioning, source locations, and selected 
connection fluxes (ipij) calculated by the MCT model as described in §3.1, for the first and 
last time. Only connections related to the seven separators discussed at length below are 
shown. In this figure, the change in the source locations between the first time and the last 
time can be compared. For instance, it can be seen that P3 and N6 are approaching with a 
rotation during ~ 24 hours indicating that this region is generating significant currents (see 
details in §3.4). 

3.3. Magnetic topology of AR 10963 

The MCT model was used to establish the magnetic topology of AR 10963 at 00:00 UT 
on July 15. The flux concentrations obtained from the MDI observations by the partitioning 
explained in the previous section were used as a set of distinct sources. Figure [6] shows seven 
of the separators for this active region, specifically those which are topologically related to 
the areas in the light curves shown in Figure [2j 

In order to track the build up of currents along the separators, it is useful to represent 
the domain fluxes (ifrij) in terms of the source fluxes (<&*.) and the separator fluxes 
where the indices i, j and k refer to sources while the index a refers to separators. Because 
all the sources are in the plane of the photosphere, there is a "mirror corona" below the 
photosphere, which is a reflection of the corona. Thus for each separator, there is a mirror 
separator, and the two form a closed curve. To express the domain fluxes in terms of the 
source and separator fluxes, it is convenient to first determine which domains pass through 
the closed curves formed by each separator and its mirror. 

To illustrate how to determine which domains pass through such a closed curve, consider 
separator 5 (hereafter \l/ 5 ) which starts and ends at the nulls Bll and A3, respectively, and 
consists of the intersection of the separatrix surfaces from these two nulls. The separatrix 
surface associated with null A3 intersects the photosphere along the fan trace from A3 to 
source P4, then along the spines of null B8 to source P2, then along the spines of null Bll 
to source P3, and returns to null A3 along the fan trace from P3. This surface encloses only 
source N6, thus all field lines within this separatrix surface must terminate on N6. Similarly, 
the separatrix surface associated with null Bll only encloses source P3, thus all field lines 
within this separatrix surface must start on source P3. Since ^5 is the intersection of these 
separatrix surfaces, the domain P3-N6 passes through the closed curve consisting of ^5 and 
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its mirror image below the photosphere. 

Alternatively, one can consider the photospheric footprint of domain P3-N6. Since all 
the sources lie in the plane of the photosphere, any field line originated in this plane must, 
by symmetry, remain in this plane. One can see in Figure [6] that any field line initiated at 
source P3 in the photosphere must pass under in order to terminate on source N6, and 
hence the domain P3-N6 passes through the closed curve of ^5 and its mirror. 

As a somewhat more complicated example, consider ^3, whose endpoints are on nulls B6 
and A2. The separatrix surface associated with B6 encloses only source P5, but the separatrix 
surface associated with A2 encloses both source N3 and source N5 (Note that source N5 is 
also enclosed by the separatrix surface associated with null Al, and this separatrix surface 
is, in turn, enclosed by the separatrix surface from null A2, but the important point is that 
source N5 is enclosed by the separatrix surface from A2). Thus, domains P5-N5 and P5-N3 
both pass through the closed curve consisting of ^3 and its mirror. Once again, one can 
instead consider the photospheric footprint of these domains. Any field lines initiated in the 
photosphere from P5 must pass under ^3 to reach source N5, and similarly for field lines 
which terminate on source N3. The domains enclosed by the separators and the nulls where 
the separators start and end are shown in Table [H 

In addition to the domains associated with each separator, each null whose fan is un- 
broken (that is, each null for which all field lines i n the fan surface end on the same source 



and thus is not the start /end of a separator; see iLongcope &: Klapperl (120021 ) ) implies the 
existence of a domain with flux equal to the flux of the spine source enclosed by this domain. 
For example, null BIO has an unbroken fan, and so the flux in domain P6-N1 equals the flux 
of source P6. This constrains those domain fluxes not associated with any separator. 

Knowledge of the domains passing through the closed curves of the separators and their 
mirrors, combined with knowledge of the nulls with unbroken fans is sufficient to express 
the domain fluxes in terms of the separator and source fluxes. The flux in each domain not 

n 

enclosed by any separator is determined from conservation of the source fluxes, $j = ?/>y, 

3=1 

where is the source flux determined from the observed magnetic field obtained by the 
partitioning process, ipy is the domain flux between the i th and j th sources, and n is the total 
number of sources of the opposite polarity to source i. The result, shown in Table HJ provides 
information about which domains can be affected by each separator. This information on 
the relationship between domain fluxes and separator fluxes is used to interpret our results 
in S4. 



3.4. Minimum Current Corona model 



In the previous section, we established the magnetic topology of AR 10963 using the 
MCT model which as sumes that the coronal magnetic field is a potential (current-free) 
field. iLongcopd ( 119961 ) developed the MCC model to incorporate source evolution, current 
generation along separators, and energy build-up in the context of the MCT model. We 
utilize the MCC model to quantify the mag netic free energy of A R 10963 in §3.5. In this 
section, we describe briefly the MCC model (ILongcopd Il996l 120011 ) . 



As a constraint in the MCC model, ILongcopd ( 120011 ) has proposed a flux-constrained 
equilibrium (FCE). This equilibrium is found by minimizing the magnetic energy subject 
to the constraint that, in the absence of emergence and submergence of flux through the 
boundary, the domain matrix does not change. This constraint is equivalent to prohibit- 
ing reconnection in the corona. The minimum energy magnetic field occurs when current 
flows only along separator field lines, such that each domain is current-free (potential) but 
is encircled by current ribbons that cancel the non-potential flux introduced by footpoint 
motions. 

In a three dimensional magnetic field, the currents occur at poin t s topologically equiv- 
alent to separators in the potential field ( iLongcope fc Cowleyl Il996l ; lLongcope fc Klapper 
20021 ). Therefore, the current in the minimum energy state can be found using the flux 
through a separator. In the minimu m energy state, the flux (^a) through a separator (a) 
can be expressed as (ILongcopd 119961 ) 



where I a is the current flowing along a separator a. The vacuum term (^a ) is the sum 
of the domain fluxes enclosed by the separator in the potential field (see §3.3). The second 
term on the right hand side, (^a^) is a self-flux from the current along the separator. The 
third term is the flux from other separator currents. 

At the first time (to), we assume that there is no current. Then, from Equation [T] we 

get 

y a (t ) = ¥:\t ). (2) 

Using the constraint on the domain flux that there is no connectivity change, Equation [1] 
becomes, at t and later time t%, 

= ^o) = <M*o = ^Hh) + ^\t x ) + J2 (3) 
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Therefore, the flux difference in the potential field (A^^ = (t\) - (to)) should be 
canceled by the sum of the second and the third terms in the right side of Equation [3] because 
the A^o- is zero with no connectivity changes. 

In order to determine the current along the separator field line, the MCC model requires 
a relations hip between the self- flux and current. The self- flux produced by current along the 
separator (lLongcopdll996l ; lLongcope fc Magard 120041 ) is given by, 



I a La- 



in 



\U 



(4) 



where e is the base of the natural logarithm, I a is the current along separator, L a is the 
length of the separator in potential field, and I* is a characteristic current proportional to 
the average perpendicular magnetic shear along separator (see Longcope & Magara 2004 for 
details). 

The third term on the right hand side in Equation [TJ the flux from other separator cur- 
rents , is expressed with a mutual inductance matrix (M^/. lLong cope!ll996l : lLongcope fc Beveridgd 



20071 ). 



E *?-' = E M °°'~< (5) 

assuming that two separators are far enough apart, so the change in the flux enclosed by 
separator a' will not be greatly affected by the internal distribution of the current on a. 
Therefore, the current along the separator field line is calculated by 



I<jL a 



In 



la' 



(6) 



The above equation is used to find the current I a , and thus the energy change for each 
separator, AWmcc, can then be calculated from the following 



AW, 



MCC 



VT,« 



2c 2 



^2 M ™' 1 °'- 



(7) 



From the above equations, we evaluate the current along each separator field line which 
arises due to footpoint motions, and the inferred excess energy over the potential field, 
referred to as the free energy. This f r ee energy is available t o be released in the form of 



localized heating (e.g., lLongcope] l200ll ; IPriest et al.ll2000l . 120051 ). 
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3.5. Magnetic energy of AR 10963 

We utilize the MCC model to investigate the magnetic energy changes of AR 10963. The 
evaluated magnetic energy changes are compared with the coronal loop brightness evolution 
observed by XRT and TRACE in §4. 

The MCC model infers the current and the free energy by using the changes of the 
potential field flux through the separators as described in the previous section. However, 
the magnetic topology evolves with ti me, so separators that are present at one time are not 



necessarily present at the next (e.g.. Ides Jardins et al.ll2009l ). Therefore one cannot follow 



the same separator for the entire time series, and a different approach is needed to determine 
the change in flux enclosed by each separator. 

The fundamental assumption of the MCC model is that there is no magnetic reconnec- 
tion in the corona, thus the flux in each domain must remain constant. As described in §3.3 
and shown in Table [TJ we have related the flux enclosed by each separator to the flux of the 
set of domains enclosed by the separator. Therefore, the flux through separators used in the 
MCC model has been replaced by the appropriate combination of domain fluxes, smoothed 
with a boxcar function with a width of 3 time steps at each time to reduce the noise due to 
partitioning. 

We evaluate the magnetic energy changes assuming that the magnetic field at the first 
time, July 14 01:39 UT, is a potential field with no current. Then, the domain fluxes are 
used to determine the changes in separator flux used in Equation [6] to derive the currents, 
and Equation [7J is used to determine the free energy. In the initial stages, the model does 
not represent the coronal currents well, since the system is started in a potential state, but it 
should improve as it evolves away from the potential state. Since the energy evolves slowly 
over several hours in this active region, the choice of different initial time within a few hours 
will not affect the results of our energy calculation. 

In Figure [7J we show the free energy resulting from current along the seven separators 
shown in the topology map in Figure [6j Separators \&i and \&2 show a steady increase in the 
free energy after ~15:00 UT. Separators ^3 and \&7 show no significant changes in energy 
from the initial potential state; thus those curves stay near the x-axis in Figure As seen 
in Figure ^3 and ^7 are nearly co-located and enclose many of the same domains, so their 
free energies behave similarly. Separator \l/ 4 shows an increase in energy until ~13:00 UT, 
and then a decrease until the last time. Separator ty§ shows an increase in energy from 
~18:00 UT and then a decrease. Separator shows a larger increase from ~18:00 UT. 

The mutual inductance term used in this analysis does not account for the internal 
distribution of current along each separator (see §3.3). Separators ^3 and ^7 are located 



-14- 



close together as seen in Figure [6], but the free energies are very small. Therefore, we assume 
that the mutual inductance is not important for the overall energetics of the system. 



4. Magnetic energy buildup and coronal loop brightness evolution 

We have evaluated the free energy build-up along separators in AR 10963 over a 24 hour 
time period using the MCC model derived from a series of 15 times of MDI observations. 
This free energy can be released through localized heating and rapid magnetic reconnection 



( )Longcopell200ll ) . We have also introduced the brightness evolution of coronal loops observed 
by XRT and TRACE. The free energy of AR 10963 is now compared to the observations in 
the Soft X-rays and EUV to determine if the free energy is related to the coronal loop 
brightening. 

In Figure HJ we show the light curves for four locations selected from the coronal loop 
observations of Hinode/XBT and TRACE. Box 1 shows a brightening in XRT shortly followed 
by a brightening in TRACE. Boxes 2, 3 show many transient brightenings in both telescopes. 
Box 4 shows steady emission in both XRT and TRACE. We compare these light curves to 
the free energy buildup for ~ 24 hours as obtained by the MCC model. 

We compare the locations of separators with the coronal loops of AR 10963 in Figure EJ 
where the magnetic topology of AR 10963 is superimposed on the XRT and TRACE images 
at 23:54 UT on July 14 viewed on the image plane, and the separators and light-curve box 
locations are also shown. In this section, we first investigate which domains lie along the line 
of sight at various coronal heights for each of the four light curve boxes. Second, we discuss 
the coronal loop brightness evolution in comparison with the domains and separators. 



4.1. Domains along the line of sight 

The optically thin X-ray and EUV emission seen by XRT and TRA CE is due to contri- 
butions from plasma along the line of sight of the observations. The MCT model is useful 
for determining which domains lie along the line of sight through the corona, at various 
heights above the photosphere. We determine which domains intersect the line of sight at 
the locations of the four light curve boxes in order to map which domains contribute to the 
observed emission. 

The left panel of Figure M shows the footprint of domains on the photosphere, along 
with the four light curve boxes (10"xl0"). Each footprint is bounded by spines and fan 
traces, shown as black thin and dashed lines, respectively. Each color represents a different 
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domain which is shown on the right side of the Figure with the same color. The right panel 
of Figure M shows potential field lines with the same colors on the left panel. The potential 
field lines are shown on the XRT image as in Figure [SJ 

Figure [10] shows the domains as a function of height along the line of sight from each 
of the four boxes where the light curves were calculated. Each box has four points on the 
horizontal axis that represent the four corners (SE, NE, SW, NW in that order) of the 
10" x 10" boxes used to calculate the light curves. The flux domains are evaluated every 
0.1 Mm from the photosphere. 

Some domains in the left panel of Figure [10] are not readily evident in the left panel of 
Figure [9] The P2-N2 domain is not seen in the left panel of Figure [9] because the domain 



does not have a photospheric footprint, i.e., it is a coronal domain ( IBeveridge fc Longcope 



20051 ). The P5-N2 domain is located between domains P5-N3 and P5-N1. It is very thin 
and can be seen only near (280", —150") in Figure [9] and near 50 Mm of Box 1 in the bottom 
panel in Figure [10] 

We show the domains for each light curve box in the second column in Table The 
domains are listed in order of coronal height from the photosphere. For instance, the light 
curve in Box 1 receives contributions from P5-N3, P5-N1 and P5-N2 below P2-N1, although 
P5-N2 occupies only a very thin layer near 50 Mm. 

To determine which domains make the main contribution to each light curve, we need to 
know the height range of the loops contributing to the light curves. We estimate the height 
of a loop by assuming that the loop is a semicircle lying in a plane perpendicular to the 
solar surface, with diameter equal to the distance between the two ends of the loop seen by 
XRT and TRACE (see §2.2 and Appendix). Using this method, we find that the loop that 
clearly contributes to the light curve in Box 1 has a height of ~50 Mm. This number takes 
into account the projection due to the location of the AR 10963 ~ 250" solar west of central 
meridian. We estimate the heights of loops that contribute to the light curves in Box 2 
- 4 (see Figure [11] and Appendix) and find that the estimated heights are all lower than 
50 Mm (see Table HJ). Therefore we assume that the heights of coronal loops that produce 
the enhancements in the light curves are lower than ~50Mm. Loops in the core of active 
regions tend to be relatively short, so this is a reasonable assumption. Using this assumption, 
we only include domains that are lower than ~50 Mm in the analysis. 

In the high corona above ~50Mm, the domains P4-N1, P7-N1, and Pl-Nl occupy the 
corona for all boxes. In addition, the domain P2-N1 extends over Boxes 1, 2, &3, but not 
Box 4. Loops contained in these high domains would mostly be very long, as can be seen in 
the potential field extrapolation in the right panel of Figure [9] Assuming these loops follow 
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RTV scaling laws ( IRosner et al.lll978f ). the density in them would be quite low because of 
their long length. Thus, we assume that the domains in the corona above the estimated 
height of the observed loops do not significantly affect the enhancements in light curves. 
The domains below ~50 Mm are represented with bold letters for each light curve box in 
Table H 



4.2. Relationships among the light curves, domains, and separators 



In this section, we investigate which separators can affect those domains that intersect 
the lines of sight for the light curves of the XRT and TRACE emission. This result is key to 
investigating if the free energy calculated by the MCC model can explain the coronal loop 
brightness evolution. 

Table [5] shows the separators that can affect each domain (see §3.3). We show the 
relevant separators next to each domain in Figure [TDJ Table [3] summarizes the relationships 
between the light curve boxes, domains, and separators. The separators that can affect the 
domain flux are shown in the third column in Table 3. As we described in the previous 
section, we consider the domains below ~50 Mm. 

The light curves of Box 1 show a transient brightening in the TRACE observation 
after a decrease of intensity in the XRT observation, implying that there is a cooling of the 
coronal loop from X-ray temperatures ( > 2 MK) through EUV temperatures (~ 1 MK, see 
Figure |2]). Observations of a coronal loop brightening in EUV following a loop brightening in 
X-rays have been examined previously. Many authors have studied the correlation between 
X-ray and EUV brightenings and the implications of this correlation for the loop heatin g 
mechanism (IWinebarger fc Warrenll2005l ; lUgarte-Urra et al.ll2009l . 120061 ; IWarren et al.ll2009l ) . 
Winebarger fc Warrenl ( 120051 ) study several loops that are relatively seclude and uniquely 
shaped so that they are able to follow the same loop as it cools through the different filters. 
They conclude that the EUV emission can be produced by a loop that emits in the X-rays 
several hours before, although the delay time is difficult to reproduce through hydrodynamic 
modeling. 

We have evaluated an approximate cooling time of the loop including the light curve 
Box 1 in Figure [TJ Assuming that plasma is at a temperature of 10 MK when it peaks in 
XRT, and is at 1 MK when it peaks in TRACE, and the loop is semicircular with a radius 
of 50", width of 10" and has a uniform cross section. With these assumptions, the density 
at the time of peak intensity in TRACE is about 2 x 10 9 cm -3 , and the loop has a cooling 
time on the order of about an hour to go from 10 MK to 1 MK. The temperature and 
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density e volution are calcula ted using an analytic solution to the ID hydrodynamic energy 



equation (ICargill et al.lll995l ). and the assumptions are not to be believed in detail. However 
this calculation is a reasonable zeroth order estimate for the cooling time. Comparing the 
estimated cooling time with the light curve of Box 1 (Figure [2]) that shows about 2 hours 
between peak in XRT and TRACE, it is likely that the loop is cooling from XRT to TRACE 
temperatures. Thus we conclude that the emission in Box 1 is probably due to a cooling 
loop. 

The light curves in Box 1 are integrated over domains P5-N3, P5-N2, and P5-N1 
below ~50Mm. These domains are related to the separators ^1,^3, and ^7. The latter 
two separators, \&3 and ^7, have relatively little free energy with very little change as shown 
in Figure Although the domain P5-N3 associated with ^1 is along the line of sight of 
Box 1, the geometry of this separator indicates that it is probably not a key component of 
the emission in this area. In the right panel of Figure [9j the loop that is included in Box 1 is 
located in a region dominated by the magnetic field lines of domains P5-N3 and P5-N1. In 
addition, the height of ^3 and ^7 are both 50 Mm, and their lengths are both about 160 Mm. 
In contrast, is much lower and shorter, with a height of 27 Mm and a length of 89 Mm. 
A few studies have found that longer coronal loops lead to longer delay times between the 



appearance of the X-ray loop and the appearance of the EUV loop (jUgarte-Urra et al.ll2006 



Aschwanden et al.ll2003l ). The delay time between the XRT and TRACE emission peaks in 
Box 1 is relatively long, thus the emission is likely to have come from longer loops that are 
located higher in the corona. The separators ^3 and ^7 are more likely associated with long 
loops than so the emission in Box 1 is probably related to separators ^3 and ^7. 

Therefore, our analysis shows that the emission variation in the area of Box 1 as summa- 
rized by its light curves most likely comes from the domains associated with the separators 
^3 and ^7, which have a small free energy. Thus it is possible that the loop is in a nearly 
potential environment, so that not much energy release can occur. This small free energy 
may be the reason why there are not many repeated transient events in Box 1, as there are 
in Boxes 2 and 3. 

Boxes 2 and 3 are related to domains P2-N3, P2-N2, P3-N3, and P2-N1. These 
domains are related to separators ^1, ^3, ^4, ^6; an d ^7. Boxes 2 and 3 display several 
transient brightenings in Figure |2j Thus it is likely that the variable emission in these areas 
comes from the domains associated with a high number of separators, including, but not 
limited to, the separators for Box 1. Additionally the free energy associated with 
and \l/ 6 increases over the time period of the observations, indicating that the fields associated 
with these separators are stressed and non-potential. This configuration could provide an 
environment for continuous impulsive events in Boxes 2 and 3 through multiple reconnection 
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processes occurring at the separators. 



Recent studies have shown that continu ous impulsive heating events are common i n 
loops located in the core of an active region fjUgarte-Urra et all 120091 ; IWarren et al.l 120071 ). 
Warren et all ( 120071 ) presented bi-directional flows as evidence of magnetic reconnection dur- 
ing thj^evojj^rtion^^ loop s in the XRT and the EIS observations. The loops observed 
by IWinebarger fc Warren! (120051 ) shift positions as they cool, which may be evidence for a 
previous reconnection event. Therefore, our conclusion that areas of active regions that show 
impulsive emission are probably related to reconnection events along separators is broadly 
consistent with previous observations. 

Transient brightenings in this active region during the observation period have different 
characteristics. The light curves of Box 3 show transient brightening in both XRT and 
TRACE. However, the light curves for Box 2 show transient brightenings in the XRT, but 
the TRACE light curve shows only very small changes. This situation may indicate that 
loops in Box 2 are being heated more frequently and have insufficient time to cool into the 
TRACE passband before being reheated. This kind of difference can not be addressed by 
our analysis since we only locate the possible sites for reconnection, and do not determine 
the frequency with which reconnection takes place. 

The light curves of Box 4, which shows steady emission in both XRT and TRACE, inter- 
sect domain P6-N1 below ~50 Mm. This domain does not have any associated separators. 
"Steady heating" has been suggested as the heating m echanism for loops that have TRACE 



moss at their footpoints, as the loop in Box 4 does (lAntiochos et al.l l2003t I Warren et al. 



20081 ). In this analysis, the emission of Box 4 comes from a domain that is not associated 
with any separators, indicating that larger-scale reconnection events are un likely to be con 



tribut ing to the heating of this plasma. Thus the steady heating proposed by lAntiochos et al. 



( 120031 ) is not likely to be due to reconnection events at separators shown in the larger-scale 
magnetic field configuration. Rather, the source of the heating that results in the emission 



observed in B ox 4 may be due to many smal 
tangling (e.g. 
MHD waves 



Antiochos et al. 



e.g. 



2003 



Kumar et al 



2006 



reconnection events or nanoflares from local field 



Parker 1983b 



a 



Ofman fc Wang 



988 



2008 



Klimchukl 120061). or altern atively, 



Antolin fc Shibatall2010h 



Several studies for the heating mechanism of coronal loops have difficulty explaining 
the brightness of observations in bot h EUV and X-ray emission and the delay between th e 
appearance of X-ray and EUV loops (IWarren fc Winebargerll2006l ; lUgarte-Urra et al.ll2006l ). 
However, recent studies have pointed out the import a nce of magnetic pro perties to explain 
the coronal loop evolutions (jUgarte-Urra et al.l 120091 ; IWarren et al.ll2009l ). Our result also 
concludes that the magnetic topology plays a significant role for the study of the heating 
mechanism of non-flaring coronal loops. 
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5. Conclusion 

In this analysis, we examined Hinode/XKT and TRACE observations of particularly 
well-resolved coronal loops in an active region that show different characteristics in their 
emission properties. In order to understand these differences, we compare the coronal loop 
evolution represented by four light curves from XRT and TRACE observations to the free 
energy calculated by the Minimum Current Corona (MCC) model. We compare the free 
energy gained from the motions of the magnetic footpoints with emission patterns that 
include steady emission and transient emission. This analysis marks the first time that the 
MCC model has been applied and tested as a full time series analysis for a quiet active region 
to investigate coronal loop evolution in plasma of different temperatures. 

This work shows that the magnetic topology of an active region can provide crucial 
information for understanding the evolution of coronal loop brightness. One of the major 
questions in coronal loop evolution studies is the difference in the heating mechanism for 
steady emission versus transient emission. We find that regions which have very dynamic 
emission are related to a high number of separators that have enhanced free energy, indicating 
that reconnections along these separators may play a role in repeated transient heating 
events. In contrast, the light curve in Box 1 that shows a single transient brightening in the 
XRT observations followed by a brightening in TRACE is related to separators with little 
free energy. This relationship indicates that the fields in the vicinity of Box 1 are nearly 
potential to begin with. Thus there is not a lot of excess energy available for release through 
repeated reconnection events in this region, and repeated transient events are not observed 
in the light curve. 

Finally, the steady emission in both XRT and TRACE comes from the domains that are 
not associated with any separator, indicating that reconnections involving the larger-scale 
field configuration are probably not associated with steady heating, and implicating either 
nanoflares or MHD waves as a possible heating source for these loops. We have found a 
topological explanation for the two classes of coronal loops, those that emit steadily and 
those that show transient emission properties. 

The Solar Dynamics Observatory (SDO) launched on 2010 Feb. 11, will undoubtedly 
shed new light on the problem of active region loop heating. The Helioseismic and Mag- 
netic Imager (HMI) onboard SDO will provide a high cadence of vector magnetogram data. 
These observations will provide detailed temporal information of the photospheric bound- 
ary evolution with which to place in context the observations by the Atmospheric Imaging 
Assembly (AIA) and its EUV observations in various coronal temperatures. The data from 
this instrument suite will complement the analysis available, as demonstrated here, from the 
Magnetic Charge Topology and Minimum Coronal Current models for further insights into 
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the mechanics of active region heating. 
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A. Appendix 

Figure [TTJ shows several loops and their ends that produce enhancements in light curves. 
Using these loops, we estimate the height of the loops that contribute to the light curves 
(see §4.1). 

Panel a) shows structures (XL4) that contribute to the light curves in Box 4. The region 
that contributes to the light curves in Box 4 is a diffuse area and does not show as a single 
loop, and there are no changes in X-ray and EUV emission. The location of the structure 
for the Box 4 is not represented in panel b) because the loop is not seen on the TRACE 
observations. An estimated height of loops using the ends of loops (F4) is represented in 
Table 1 (see also §4.1). 

The transient brightenings in the light curve Box 2 and Box 3 are produced by several 
different loops at different times. Therefore, in panels c)-h), we show the loops that produce 
enhancements in the light curves of Box 2 and Box 3 at three different times when the 
transient brightenings are observed by XRT and/or TRACE. In panel c), each loop appearing 
in Box 2 (XL2) and Box 3 (XL3') that shows enhancements in the light curve in the XRT 
observation at around 20:00 UT is represented. The ends of loop seen by XRT for Box 2 is 
marked by F2 in panel c). The loop (TL3') that shows enhancements in the TRACE light 
curve in Box 3 is a small loop shown in panel d). In panels e) and f), we show the loops 
that correspond to enhancements at around 22:00 UT. The ends of loop seen by XRT and 
TRACE for Box 3 are marked with F3 and TF3 in panels e) and f), respectively. In panel 
e), the loop (XL2') that contributes to the light curve Box 2 at that time is also shown. 
In panel g), the loop (XL2") that shows the enhancement in the light curve for Box 2 at 
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around 01:00 UT is shown. At this time, there is no distinct loop in Box 3 in the XRT 
observation, although there is a small enhancement in the light curve. In panel h), no loops 
are represented because there is no enhancement in the light curves and brightening in the 
TRACE observations at this time. The loops appearing in Box 2 and Box 3 at the last three 
times in Figure ITT1 look like different loops (see animation movie), though their locations are 
similar. 

We represent the height of loops, LI, XL2, XL3, TL3, and XL4, in TablelH The distance 
between the ends of loops, XL3', TL3', XL2', and XL2", are similar or smaller than the loops 
in Table HI Therefore, the heights of these loops are lower than the heights of the loops in 
Table H 
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Table 1. Separator properties 



Separators 


Nulls 


Length (Mm) 


Height (Mm) 


Domain 


*l 


A1/B6 


89 


27 


(P5-N5) 


*2 


A2/B11 


55 


8 


-(P3-N3)* 


*3 


A2/B6 


153 


50 


(P5-N5)+(P5-N3) 


*4 


A3/B8 


104 


22 


-(P4-N6)* 


*5 


A3/B11 


44 


11 


(P3-N6) 


*6 


A4/B11 


61 


11 


(P3-N6)+(P3-N1) 


* 7 


A4/B6 


160 


50 


(P5-N5)+(P5-N3)+(P5-N2) 



* ' — ' represents an opposite direction in flux flows between flux domain and separator 



Table 2. Connectivity matrix in terms of source and separator fluxes 





Poo 


PI 


P2 


P3 


P4 


P5 


P6 


P7 


Nl 




$P1 


* 


*6 - *5 


*P4 + *4 


<I> P5 _ y 7 


$P6 


$ P7 


N2 








*JV2 - *P3 + *6 - *2 - *7 + *3 


<S>P3 ~ *6 + *2 





* 7 - * 3 








N3 








$N3 +*2 - *3 + *1 


-*2 





*3 - *i 








N5 
















*1 








N6 








*JV6 - *5 + *4 


*5 


-*4 












Note. — Poo is a positive source "at infinity" introduced to balance the flux from active region. The flux of P c 
is about 10 % of total unsigned flux. 

*$JV1 - $ Poc - $P1 - <I>P6 - *P7 - *6 + *5 - *P4 - *4 - *P5 + *7 



Table 3. Domains and separators that influence light curve boxes 



Boxes 




Flux domains 


Separators 


Box 1 


P5- 


-(N3, N2, Nl), P2-N1, P4-N1, P7-N1, Pl-Nl 


*1, *3, ^7 


Box 2 


P5- 


-N3, P2-(N3, N2, Nl), P4-N1, P7-N1, Pl-Nl 


*1, *2, *3, *4, ^6: *7 


Box 3 


P3- 


-N3, P2-(N3, N2, Nl), P4-N1, P7-N1, Pl-Nl 


*1,*2,*3,*4,*6,*7 


Box 4 


P6- 


-Nl, P4-N1, P7-N1, Pl-Nl 





Note. — Bold represents the domain that used to relate separators. 



-26 - 



Table 4. Loops and estimated heights 



Loops 


Time 


Height (Mm) 


Note 


LI 


22:24 UT* 


50 


Figure [JJ 


XL2 


20:07 UT 


40 


Panel c) in Figure [TTJ 


XL3 


22:07 UT 


48 


Panel e) in Figure [TT] 


TL3 


22:07 UT 


14 


Panels f) in Figure [TTJ 


XL4 


22:24 UT 


46 


Panels a) in Figure [TTJ 



Note. — '*' represents a time when the XRT observation shows 
a brightening. 




Fig. 1. — The locations of four boxes for light curves (Figure |2]) on XRT (left) and TRACE 
(right) images. Smaller and larger boxes represent the areas of ~ 5" x 5" and ~ 10" x 10", 
respectively. A Loop (LI) that produces the enhancement in light curve Box 1 is represented 
by "+" symbols along the loop. The ends of this loop are represented by Fl. This figure is 
also available as a mpg animation. The movie shows the XRT and TRACE observations for 
~ 6 hours corresponding to the light curves in Figure [2J 
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22:00 
Time 
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Fig. 2. — Light curves for the four boxes in Figured) Black and blue lines (solid: ~ 5" x 5", 
dotted: 10" x 10") represent XRT and TRACE, respectively. Left and right axes repre- 
sent the DN s _1 arcsec^ 2 for XRT and TRACE, respectively. Left top (Box 1): A bright- 
ening in TRACE following a brightening in XRT, indicating a single cooling event. Right 
top (Box 2): Transient brightenings in XRT and only very small changes in TRACE. Left bot- 
tom (Box 3): Transient brightenings in both XRT and TRACE. Right bottom (Box 4): Steady 
emission in both XRT and TRACE. 
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Fig. 3. — Top: Partitioning at 00:00 UT on July 15. Pole locations of positive (P: +) 
and negative (N: x) sources are represented. Contour represents the partitioning. The 
background MDI images show only the smoothed field above the threshold used in the 
partitioning algorithm. Bottom: Colored lines represent domain fluxes among all sources. 
The contours surrounding individual pixels which are included or excluded from a partition 
are rendered as small diamonds. 
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Fig. 4. — Total unsigned source flux for all sources obtained by the time-series partitioning. 
The integration times of the MDI observations are represented with uncertainties in the total 
unsigned flux: 300 sec (x) and 30 sec (+). Dotted line represents the total averaged flux 
used in this analysis excluding the source flux of N4 ($tv4 — 0, see details in text). Averaged 
uncertainty for 15 times in the total unsigned flux is 5.0 x 10 19 Mx. 
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Fig. 5. — Seven domains related to the separators (Table [I]) are shown on the top and the 
bottom panels at 01:39 UT on July 14 and 00:00 UT on July 15, respectively. Colored lines 
from each positive (+) and negative (x) source represent the flux in each connection. The 
scale of color is the same in both top and bottom panels, but this scale is different from the 
right panel of Figure |3j Footpoint motions and the changes of domain fluxes over the 24 
hour period can be seen by comparing the top and the bottom panels. 
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Fig. 6. — Topology of AR 10963 at 00:00 UT on July 15. Right panel shows an enlarged 
version of the middle of left panel. We show the topology on the tangent plane to explain 
the domains enclosed by separators in §3.3. Nulls (B-type: red y and A-type: blue A) 
are located between positive and negative sources, respectively. Black dashed and black 
thin solid lines are fan traces and spines that show where separatrix surfaces intersect the 
photosphere. Colored lines are the separators: the intersection of separatrix surfaces. ^3 
(A2/B6) is represented by a thinner line than other separators to be visible with $7 (A4/B6) 
because the locations of these two separators are close except between nulls A2 and A4. 
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Fig. 7. — Energy buildup from 01:39 UT on July 14 (current=0 at this time) for each 
separator. Colors are the same as in Figure [6j Black thick line represents the start time of 
the light curves in Figure [2j ^ and ^2'- A steady increase after ~15:00UT. \l/ 3 and ^7: no 
significant changes in energy from the initial potential state. ^4: an increase until ~13:00 UT 
and then a decrease. ^5: an increase from ~18:00UT and then a decrease, an increase 
from ~18:00UT. 
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Fig. 8. — Locations of separators on XRT and TRACE observations. The symbols are 
the same as in Figure [6j Fan traces (white dashed lines) , spines (white thin solid lines) , 
and separators (colored lines) are placed on the view from the Earth (image plane). The 
separator colors are the same as in Figures |H1 H The locations of separators are slightly 
shifted comparing with the locations in Figure |6] due to the difference between the image 
plane and the tangent plane. Light curve boxes are also represented as four boxes. 
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Fig. 9. — Left: Domain footprints on the image plane. Symbols are the same as in Fig- 
ures El [HI The domain footprints are divided by the spines (black solid lines) and fan traces 
(black dashed lines) on the photosphere. Domains present along the lines of sight for the 
light curve boxes are shown in color. The positive and negative sources of the domains are 
labeled on the right side of the figure. Right: Potential field lines of each domain with the 
same colors as in the left panel. White boxes are the light curve boxes. The loops in Figured] 
are also represented with "+". 
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Fig. 10. — Domains along the line of sight through the corona at the locations of the light 
curve boxes. The top and bottom panels show the domains over a large height range and a 
small height range, respectively. The domains are shown with the same colors as in Figure EO 
Each box has four points in the horizontal direction that represent the four corners (SE, NE, 
SW, NW in order) of ~ 10" x 10" square of the boxes. The separators that could affect 
the domain fluxes are represented on the right of the domains with the same colors as in 
Figure El 
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Fig. 11. — Loops that produce the enhancements in light curves are represented by "+" 
symbols along the loop, a) XRT image for loops that contribute to the light curve in Box 4 
b) Corresponding TRACE image. c)-h) Loops that produce the enhancements of the light 
curves Box 2 and Box 3 at three different times in Figure [2j The ends of loops are shown 
as a symbol "x" with letters, F2, F3, F4, and TF3 for several loops, which are applied to 
estimate the heights of the loops (see Table H|). Light curve boxes are also shown as four 
boxes. 



